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ABSTRACT 

A satisfactory theory for the formation of globular star clusters (GCs) has long been elusive, perhaps 
because their true progenitors had not yet been guessed. In this paper I propose a causal relationship 
between the strongly decreasing densities of Lyman a clouds at high-redshift and the formation of 
globular clusters — namely, that GCs were created by the accretion of Lyman a clouds onto protogalaxies. 
I describe a scenario which involves an inherently stable and orderly cycling of compression and cooling 
in the central cores of these clouds during the extended period of dissipation in the outer regions of gas- 
rich protogalaxies, culminating in a burst of efficient star formation. I demonstrate that the comoving 
density of GCs is comparable to that of Lyman a clouds at high-redshift, that the energetic requirements 
for compression to core GC densities can be met, and that the time-scale for cooling is within obvious 
limits imposed by dynamical stability. 

This dissipative process requires there to be a large column of dissipated gas about the attractor in 
order to form GCs. In adition, the energy requirements of compression to GC central densities requires 
attractor masses greater than that capable of sustaining circular velocities of ~ 40kms~ 1 . If this scenario 
is supported by numerical simulations, then by implication, the GCs were formed at modest redshifts of 
z ~ 1 — 3. This knowledge could help to break the troublesome degeneracy between look-back time and 
redshift. This model is consistent with a picture of hierarchical galaxy growth over time scales of many 
billions of years. 

Subject headings: galaxies: evolution — galaxies: formation — globular clusters: general — stars: 
population II 



1. INTRODUCTION 

Globular clusters (GCs) are at once objects of interest 
to students of the Galaxy and to cosmologists — tothe 
former, as probes of galaxy formation and as dynamically 
pure objects of curiosity, and to the latter because of their 
great age, and for what they may tell us about conditions 
in the early Universe. Attempts to explain them have 
included models based on the accretion of protogalactic 
clumps (Searle & Zinn 1978; Searle 1980), two-phase col- 
lapse (Fall & Rees 1985; Murray et al. 1993), and cloud 
collisions (Murray & Lin 1989; 1993; Harris & Pudritz 
1991; McLaughlin & Pudritz 1996). In this paper I sug- 
gest a dissipative scenario in which Lyman a clouds are 
viewed as the progenitors of the GCs. 

Studies of high-redshift Lyman a clouds (e.g., Lu, 
Wolfe, & Turnshek 1991; Bechtold 1994) show that the 
high line density of clouds at z » 3 decreases rapidly 
with time. A major constituent of the Universe at high- 
redshift, these clouds are not likely to have disappeared 
without a trace. Binney (1976) first suggested that the 
earliest star formation in a collapsing protogalaxy would 
be within transient sheets of compressed gas generated by 
the collision of infalling gas. Much of the halo may have 
been created in this way. But the formation of GCs re- 
quires a more extraordinary variant of this accretion pro- 
cess. If clouds are centrally condensed (e.g., dark matter 
(DM) held), this would help them to survive the infall 
intact, possibly leading to a super-concentrated burst of 
star formation, and hence possibly a GC. The proposed 
GC formation scenario utilizes the kinetic energy of infall 



to compress clouds to densities comparable to the central 
regions of GCs. 

Below, in referring to galaxian quantities and galacto- 
centric distances, I use upper-case letters, and for cloud 
quantities, lower-case letters. When this convention isn't 
applicable, I will use subscripts. 

2. ANECDOTAL EVIDENCE 

Zaritsky (1995) has shown that there are abundance gra- 
dient anomalies, and asymmetries in the distribution of HI 
at large galactocentric radii in 1/3 to 1/2 of nearby spi- 
ral galaxies. This was cautiously interpreted as being the 
result of the accretion of low metallicity HI clouds within 
the last few Gyr. More recently a measure of lopsidedness 
and anisotropy has been developed from the Fourier anal- 
ysis of the surface brightness distributions of disk galaxies 
(Zaritsky & Rix 1997; Rudnick & Rix 1998; Jiang & Bin- 
ney 1998) which has been found to correlate well with the 
locations of enhanced star formation and H I anomalies. 
These findings support the conclusions of Zaritsky's (1995) 
earlier work. 

A source for these apparent accretion events may have 
been found in a population of low redshift HI clouds, 
discovered in Hubble Space Telescope (HST) spectra of 
low redshift quasar and active galaxy spectra (Bahcall et 
al. 1991; Morris et al. 1991; Morris et al. 1993, Lanzetta 
et al. 1995; Stocke et al. 1995; Shull, Stocke & Penton 
1996; Tripp, Lu, & Savage 1998), and at redshifts out to 
z ss 0.6 (Chen et al. 1998). Many of these clouds appear 
to be clustered about luminous (My <, —19) field galaxies. 
Recently, Tripp et al. (1998), adding data from two QSO 
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sightlines, have summarized recent developments, showing 
that of those clouds which are found within 2/i 7 T 5 1 Mpc of a 
major galaxy, most have a projected distance of less than 
700 kpc. However, nearest-galaxy distances are seen to be 
greater than 2h^ Mpc for as many as 40% of the clouds, 
prompting Tripp et al. (1998) to refer to this population 
as "void clouds" . 

In a study of high velocity clouds (HVCs) in the Local 
Group, Blitz et al. (1998) note that the internal veloc- 
ity dispersions of remote HVCs are strongly peaked, with 
a w 20 ± 6.5 kms -1 , which argues for a generally homo- 
geneous population that is only mildly subgalactic in its 
virial temperature. Their dust-to-gas ratios are at least 
a factor of three below that of normal Galactic clouds, 
and heavy -clement abundances well below solar values. 
Wakker et al. (1998) have found that the "C" cloud, the 
closest HVC, has a metallicity, Z = 0.07 ± 0.02 Z Q , well 
below what would be expected for a Galactic fountain or 
tidal tail. I therefore assume that these apparently in- 
falling clouds are indeed representatives of the extragalac- 
tic Lyman a population. 

It therefore seems likely that the more clustered cloud 
population is a major source of the impactors hypothe- 
sized by Zaritsky (1995), and others. They apparently 
share the same physical region occupied by dwarf galaxies 
that cluster about dominant field galaxies. If clouds share 
the kinematic characteristics of dwarfs as well, then we can 
expect that clouds will cluster about their primary within 
velocities of < 400 km/s (Zaritsky et al. 1997, and refer- 
ences therein). One reason why this might be expected 
is if the dwarf galaxies, and the clouds, were both repre- 
sentative of a single spectrum of systems that inhabit the 
neighborhood of giant field galaxies. 

It is reasonable to propose that if there is a continu- 
ity between the nature of low, and high-redshift Lyman 
alpha absorption systems, then the large population of 
HI clouds at high-redshift may likewise have been loosely 
clustered about protogalaxies, and thus have been sub- 
ject to accretion to the more dominant gravitational po- 
tential in their midst. At modest rcdshifts (z <~ 0.5), 
Chen et al. (1998) has found that clouds of Lya equiv- 
alent width W w 0.3A lie at a characteristic distance of 
160/i _1 kpc (h = i^o/lOO kms^Mpc -1 ), consistent with 
more local data (Lanzetta et al. 1995; Tripp et al. 1998). 
It does not seem reasonable that these clouds are rota- 
tionally supported, for the velocity dispersion of galaxies 
within <~ 5/i _1 Mpc of the Local Group is on the order of 
only 50- eOkms" 1 (Giraud 1986; Schlegel et al. 1994), 
so that clouds of extragalactic origin with tangential ve- 
locities of ~ 200 km s -1 would seem to be highly unlikely. 
For clouds not rotationally supported, however, the infall 
time is only ~ 1/4 that of the look-back time to z <~ 0.5. 
Therefore I conclude that these clouds, and those at higher 
rcdshifts, must be being accreted to the galaxies about 
which they cluster, and further, the distribution must be 
being replenished to some significant degree since roughly 
the same distributions maintain atz« 0.6 as locally, as 
we have seen. 

At rcdshifts in the range of 2 to 3, Fernandcz-Soto et 
al. (1996) have shown that there is strong clustering within 
250 km s -1 among clouds that have associated metal lines, 
while in the range, 2.7 < z < 3.7, Chernormordik (1995) 



has found weak cloud clustering within velocity scales of 



150 km s for clouds with columns N > 10 cm 



Us- 



ing cross correlations of absorption systems in double and 
group quasars at redshifts of ~ 2, Crotts & Fang (1989) 
found that, Lyman a systems with Wo > 0.4A cluster 
on scales as large as 0.7 (1.2)/i _1 Mpc for q Q = 0.5(0.0), 
though there is an apparent significant variation with en- 
vironment. 

These findings are broadly consistent with the physical 
association of clouds with galaxies seen at low redshift, and 
we may therefore postulate that high-redshift clouds may 
accrete to their more massive neighbors on time-scales of 
~ 1 + Gyr, their dissipated remains precipitating bursts of 
star formation. 

2.1. The Density of GCs and Clouds 

If Lyman a clouds are to cause the formation of GCs at 
high redshift, then the comoving density of Lyman alpha 
clouds at high-redshift should be comparable to the den- 
sity of GCs locally. The density of GCs can be calculated 
from the luminosity function and the luminosity-weighted 
specific globular cluster frequency, Sn = N t /Lis (Harris 
& van den Bergh 1981), where N t is the number of GCs in 
the galaxy, and L15 is the luminosity in units of a galaxy 
with My = —15. The number density of GCs is given 
by the integral of the median specific GC frequency times 
the luminosity-weighted Schcchtcr (1976) function. Given 
the parameterization, X = L/L*, we find the comoving 
density of GCs is, 
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= (S N )L* 15 J X<S>(X)dX w (S N )L* 15 $*, 



(1) 



where the result follows from setting a = — 1. I make 
the following substitutions: For a local sample of galax- 
ies (central-cluster ellipticals excluded), (Sn) ~ 4 (Har- 
ris 1991), L* 15 w 100 if My w -20.0, and $* ~ 
0.014 /i~ 3 Mpc 3 (Loveday et al. 1992). Thus, 



n gc rts 5.5 h 3 Mpc 



(2) 



For clouds at high-redshift, the attributed comoving 
density of clouds of radius r c i, and line density dN/dz is, 

ncl = 17 ^r\(i+z)z ) ■ 1 assume a l° w ^ FKW cosmology 
(A = 0), for which dz/dR + z) 2 , where H is the 

Hubble constant. This yields, 



n c i 



dN/dz 



kRotKI + zY 



(3) 



where R = c/H D . A closed model doubles the expected 
density when z — 3. 

We must decide what fiducial value to substitute for 
r c i . A numerical simulation of high redshift self shielded 
primordial clouds, in which cold, isothermally distributed 
DM comprises 90% of cloud mass (Manning 1992), showed 
that the clouds responsible for most absorption lines were 
large, just short of the size at which thermal instabil- 
ity from self-shielding would occur in the central re- 
gions. These clouds were found to have H I columns of 
LogiV = 14(13) at projected radii of - 20 (45) /i -1 kpc, 
and core radii on the order of 2 kpc, and core bary- 
onic mass of <~ 1 — 4 x 10 6 M Q ). They have idealized 
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thermally broadened velocity dispersions in the range of 



17 - 22kms 



Meanwhile, Shull et al. (1996) have 



estimated the radii of local clouds to be ~ 100 kpc at 
columns of N ~ 10 13 c m - 2 , and Blitz et al. (1998) have 
found the upper limit cloud radius, constrained by the 
tidal field of the Local Group, to be ~ 25 kpc at columns 
of N ~ 10 19 cm -2 . While these figures may make the sim- 
ulated cloud sizes at high-redshift seem low, the 50-fold 
drop in the far-UV metagalactic flux at the Lyman limit 
from its value at z = 2.5 (Haardt & Madau 1996) will help 
to account for this apparent disparity in absorption cross 
section. In a study of the lensed quasar Q2345 + 007A,B, 
Foltz et al. (1984) have placed a lower limit on the char- 
acteristic diameters of clouds at (z) = 1.95 of 5 — 25 kpc 
(the uncertainty stems from the unknown redshift of the 
lens). These values simply reflect the range of possible 
distances between beams. However, a Monte Carlo anal- 
ysis of the correlated absorption systems in the spectra of 
Q1343 + 2640A, B (Dinshaw et al. 1994), which has a simi- 
lar ratio of "hits" to "misses" (h/m~ 2), found an inferred 
cloud diameter 2.5 times that of the median beam separa- 
tion. If this correction is applied to the Foltz et al. (1984) 
data, the cloud range in radii becomes 7.5 — 62.5/i _1 kpc. 
However, recent work correlating absorption systems of 
double and group quasars (e.g., Fang et al. 1996; Crotts 
& Fang 1998) find much larger values (typically fewxlOOs 
of kpc), but it is not clear that these are large individ- 
ual clouds rather than the larger-scale clustering of much 
smaller individual clouds that we have been discussing, 
which could be expected to give "cloud diameters" of from 
~ 300 — 500/i _1 kpc, depending on the equivalent widths 
probed and the method of analysis - consistent with their 
results. 

For high redshift clouds (z ~ 2 — 3) I adopt a fiducial 
radius of 25 /i _1 kpc for individual H I clouds at columns of 
Log N ~ 14. A line density of dN/dz w 100 at a redshift of 
z <~ 3 is assumed, based on data normalized to an equiv- 
alent width limit of W a = 0.24A (Weymann et al. 1998), 
which, for T w 2 x 10 4 K, is quite close to a column density, 
N = 10 14 cm~ 2 . Substituting these values into Eq. 3, we 
find, 

n cl (z = 3) w 4.2h 3 Mpc~ 3 . (4) 

For = 1, n c i = 8.4h 3 Mpc~ 3 . Comparison of Eqs. 2 
and 4 shows n gc w n c ;. If there is a spectrum of cloud 
sizes (I find the largest dominated the absorption line 
density), then since cloud number density, n c i oc r~ 2 (Eq. 
3), the total number density of clouds may be significantly 
larger, which may help the numbers to look better for open 
cosmologies. Also, if replenishment is occurring at high- 
redshift by the condensation of clouds in voids and sub- 
sequent movement toward neighboring galaxies, then the 
total number of potential impactors could be substantially 
greater than estimated above. 

3. THE PHYSICAL REQUIREMENTS FOR CLOUD SURVIVAL 

For a modest sized protogalaxy, the infall velocities of 
an intergalactic cloud can be expected to be in excess of 60 
km/sec. In order to form a cluster, therefore, some small 
fraction of the plunging cloud must survive the supersonic 
shock in a compact form. Numerical simulations of homo- 
geneous clouds subjected to interstellar shocks have found 
that clouds are destroyed on time scales of a sound cross- 



ing time (e.g., Klein, McKee & Colella 1994; Murakami & 
Ikcuchi 1993). However, numerical studies of the survival 
of intergalactic, cold dark matter-held, "mini-halo" clouds 
(Rees 1986) subjected to supersonic flows (Murakami & 
Ikcuchi 1994) have shown that the cores of clouds confined 
by a dark halo may survive extended periods of exposure 
to supersonic wind as long as the central cloud density ex- 
ceeds that of the ambient medium through which it passes. 
In the context of accretion to a galaxy, this will require 
that the central density of the cloud should increase in 
time in order that it may survive the increasing densities 
it will encounter as it falls inward. This in turn will require 
cooling. For a stable compression process, therefore, the 
cooling time scale, r c , must remain smaller than the time- 
scale for the increase in the ambient protogalaxy gas den- 
sity, Td, at each stage of the cloud's journey. The timescale 
for change of density is expressible by the equation, 



Td 



1 dp 
~p~dt 



(5) 



I presume, for the sake of simplicity, that the gas and the 
DM arc distributed in an isothermal profile with density 
law, p = IC/ (AttR 2 ), where R is the galactocentric distance, 
and K. is a unit of mass per unit length. JC will be referred 
to as the system mass-distribution constant. For an L* 
galaxy, fC = fC* w 1.14 x lO 7 M pc _1 . I also presume 
that the distribution is truncated at a distance R t . Then 
for R < R t , M(R) = ICR, and for R > R t , M(R) = tCR t . 
Manipulating Eq. 5, we find, 



Td = 



R_ 

2n~ 



(6) 



where Vi is the infall velocity. With the assumed mass 
distribution, we can calculate the potential: 



$(R < R t ) = -GK In 



1 



Rt_ 

Ro 



(7) 



where Ro is the location, presumed in this case to be <~ 2 
Mpc, of the turn-around radius, and at which the poten- 
tial is set to zero. R t is assumed to be 500 kpc for an 
L*-type object. Conservation of energy requires that the 
infall velocity be given by the equation, V{(R < Rt) < 
— 2$(i? < R t ) the inequality stemming from dissipative 
effects. Substituting this into Eq. 6, and using the value, 
JC* (above), we find that the timescale for change of den- 
sity is, 

Rkpc 



T d > 1.38 x 10 e 



160 



yr, 



(8) 



where Rk pc is the galactocentric radius in kiloparsecs. At 
a distance of 15 kpc Td ~ 1.3 x 10 7 yr. 

Dissipation is expected to begin in earnest at an esti- 
mated galactocentric distance of <~ 160/i kpc, for that is 
the projected distance at which lines of sight encounter 
EWs of - 0.3 1, or columns of - 10 14 cur 2 (Chen et 
al. 1998). It will be important that the rate of cooling 
keep pace with the time-scale for compaction. I assume 
that dissipation balances the rate of change the cloud po- 
tential energy so that the velocity remains at a constant 
SOkms- 1 . Thus we have, 



GM 



Pcl dV = n 2 kT/3 B f dV, (9) 
v Jv 
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where we will insert values for an L* galaxy, and cloud 
baryonic density, p c \ — nfimn. The case B recombi- 
nation cooling coefficient may be expressed as, (3b = 
9.17 x 10" 14 (T/20000)°- 5 cm 3 s" 1 . From this relation and 
Eq. 9, we derive, 



r= 1-127 x lQ4 / 15_ 



n 2 / 3 \Rkpc, 
The cooling time-scale of the cloud can be defined, 



A ' 



With A = n 2 kTf3 Bl we hnd 



2.313 x 10 11 
nT 1 / 2 

Combining Eqs. 10 and 12 we find, 
6.9 x 10 5 / R kpc 



t c = 



n 



2/3 



15 



yr. 



(10) 



(11) 



(12) 



(13) 



The pressures behind the head of the bow shock will re- 
sult in baryon densities p cl w 4:M 2 p ga i(R)(T gal /T cl ) dur- 
ing the comprcssonal stage, where M. is the Mach number, 
expected to range up to w 10. Assuming an L* galaxy, at 
R w 15kpc, and assuming 90% DM, then the galactic 
baryon number density is, n ga i ;> 10~ 2 . Thus, we may ex- 
pect n J> 1, and should eventually reach densities of order 
10 4 cm~ 3 . This will assure that r c is much less than r<j. 

For galactocentric radii of interest, the cooling time scale 
within the core of the cloud appears to be comfortably less 
than Td- During the final stages of accretion, another term 
must be added to Eq. 9 to represent the deceleration of 
the cloud. This will compensate for the decline of the first 
term, leading to only modest changes in the derived cool- 
ing time-scale. It must be emphasized that there should 
be a comfortable gap between Td and r c , for we don't ex- 
pect real galaxies to have entirely smooth density profiles. 

4. THE MECHANISM FOR CLOUD COMPACTION 

In order to satisfy the requirements that these events 
produce globular clusters, the cloud must be compressed 
to densities comparable to, or in excess of that found in 
globular cluster cores. Powered by the kinetic energy of 
the cloud, the combination of pressurization and radiative 
cooling is capable of compacting the cloud. In this sec- 
tion I discuss three important facets of the GC formation 
scenario. 

4.1. Dynamical Features of the Interaction 

Let us now suppose that a centrally condensed cloud 
(core ^3 — 5 kpc diameter, is drawn into the potential well 
of a gas-rich protogalaxy. When the cloud begins to en- 
counter the dissipated gaseous halo at velocities of <~Mach 
10, a bow shock is established, forming a shocked shell 
about the cloud. In the frame of the shock, galactic gas is 
quickly decelerated by the shock within this shell, impart- 
ing a momentum to the infalling cloud. Only gas of density 
comparable to that of the wind is stripped — a result of the 



small mean free path in relation to cloud size — while the 
dense core is found to survive in excess of 3 Gyr (Murakami 
& Ikeuchi 1994). Accordingly, we should expect only the 
core to survive intact, carrying <~ few x 10 6 M Q of baryons, 
as noted in §2.1, in the range of 10~ 2 — 5 x 10~ 4 of the total 
cloud baryonic mass. The shocked and stripped remain- 
der, much of which would be caught up in eddies, may 
participate in star formation in the manner suggested by 
Binney (1976), contributing to the halo population at rel- 
atively low efficiencies. The remnants of this would settle 
toward the disk. 

While the DM component is thought to play an im- 
portant role in containing the cloud in the IGM (Blitz et 
al. 1998; Shull et al. 1996), once the ram pressure, inte- 
grated over the surface of the cloud, exceeds the maximum 
force between the baryonic cloud and the DM cloud, the 
DM halo would be lost. This is expected to occur in the 
early stages of intense dissipation, at distances less than 
<~ 160/i _1 kpc (see §3). When the deceleration on the cloud 
begins to dwarf its self-gravity, the cloud will shift and 
compress to approximate an exponential atmosphere with 
a scale-height determined by the deceleration rate, which 
in turn will heat the cloud. 

Within this environment, two forces act upon the cloud 
which tend to induce vorticity, a shift of the higher den- 
sity gas within the cloud toward the shock front, displac- 
ing lower density gas toward the side, and the shear force 
due to the stripping flow inside the bow shock. While by 
itself, molecular viscosity is not sufficient to appreciably 
effect the vorticity, a turbulent viscosity will be produced. 
The shock will stand off from the cloud due to the high 
pressures immediately behind the shock, but the dense, 
cool cloud immediately "downstream" of this region will 
be subject to Rayleigh-Taylor instabilities since cloud de- 
celeration will result in denser material overlying less dense 
material. As these perturbations grow they will become 
subject to the shear flow closer to the shock, transforming 
them into Kelvin-Helmholtz (K-H) instabilities. These in- 
stabilities would continue to grow along the "fetch" of the 
interface of galactic and cloud gas, causing a turbulent vis- 
cosity and an increased transferral of angular momentum 
to the cloud. The combined effect of these forces would be 
to induce vorticity into a torroidal volume whose axis is 
aligned with the vector of the cloud's motion. During the 
stage of intense dissipation, it is better to compare it to an 
inverted convection cell than to a vortex ring, for unlike 
the vortex ring, the cloud is contained from without by the 
shock and by the deceleration, rather than from within by 
a vacuum. 

During the final stage of intense dissipation, the in- 
creased deceleration means that gas which has been 
stripped need not be irrevocably lost to the cloud. Two 
factors are responsible: 1) stripped material will be pushed 
toward the central axis of motion at ^sonic speeds because 
of the near vacuum behind the cloud, 2) the deceleration 
of the cloud will result in a relative (negative) accelera- 
tion between the stripped material and the shock front. 
The re-accretion of this gas to the cloud would reinforce 
the pattern of motion established in the early adjustment 
of the cloud. It is an important feature of this scenario 
that the stripped, and eventually re-incorporated material 
will contain metals present in the protogalaxy, which have 
been mixed with cloud gas by K-H instabilities. The well- 
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established observational constraint that the metallicity 
in almost all GCs is uniform to within 0.05 dex (Suntzcff 
1993), requires that turbulent mixing should occur well 
before star formation begins (Murray and Lin 1993). In 
this scenario, the enhancement of the metallicity, and the 
mixing, occur as a result of the natural dissipative process 
of accretion outlined above, comfortably preceding frag- 
mentation (see §4.4). 

4.2. The Stability of the Plunging Cloud 

The modeling of intergalactic clouds held by DM (see 
§2.1) indicates they may be centrally condensed. Further- 
more, centrally condensed clouds have been found to sur- 
vive shocks (Murakami & Ikeuchi 1994). Therefore, if this 
model of Lyman-aclouds is correct, they should be ex- 
pected to survive the gentle onset of ram pressure. The 
vorticity induced by the shear flow inside the shock does 
not appear to be a destabilizing influence, because the 
cloud is contained by the shock shell and by deceleration, 
as described above. However, when the cloud decelerates 
to subsonic velocities, it will evolve toward an ordinary, al- 
biet large, vortex ring. Here, the orderliness of the motion 
will be critical in determining the length of time the cell 
should survive. Studies of vortex rings at relatively high 
Reynolds numbers show that sinusoidal bending modes 
may grow, and eventually destroy the ring (Widnall & Sul- 
livan 1974). However, these modes, whose number are few 
when vorticity is widely distributed, would be damped by 
the shock which wraps around the cell during the period 
of intense dissipation, but may grow after it is no longer 
supersonic. At the very minimum, the cell will survive 
for an eddy turn-over time, given that the ordered motion 
will produce a vacuum in a ring which will cause the self- 
induced motion characteristic of vortex rings. For a cell of 
dimension 10 pc, the eddy turn-over time is, 

r e w 7.6 x 10 5 (^) yr, (14) 

where v is the peak tangential velocity of the vortex cell. 
The self-inducing motion of the vortex ring is likely to 
cause it to last many times greater than this. However, 
accurate estimates of this must await a careful numerical 
simulation, as energy and time-scale arguments are held 
hostage to uncertainties in the levels of turbulent viscosity 
during the stage of intense dissipation. 

4.3. Energetic Requirements for Cloud Compression 

If radiative cooling can keep pace with heating, we may 
assume isothermal pressurization. If the cloud core has a 
baryonic mass of m c i, an initial baryonic density of pi, and 
a final density of pj, the required work for compression is, 

W = m c ic 2 s ln{^j. (15) 

Values for W are insensitive to variations by factors of a 
few in the initial or final densities. Equating this work 
to the kinetic energy of the cloud core we find that the 
required infall velocity v i7 is given by, 



Vi = c s J2ln(p f /p t ), (16) 



where c s is the sound speed of the cloud core. For a fi- 
nal density equal to that of the median Galactic globular 
cluster core densitiy, <~ 10 3 M©pc~ 3 (Djorgovski 1993), a 
factor of compression in the range of ~ 10 6-8 is expected, 
and yields a required infall velocity, Vi/c s ~ 5.0 — 6.0. For 
a temperature of 2.0 x 10 4 K, the minimum required infall 
velocity is ~ 60 — 75 km/s, a velocity attainable even in 
dwarf systems (see e.g., Wyse & Silk 1985). 

4.4. The Transformation of Cloud to Cluster 

The final stage of the transformation is reached when 
the velocity approaches subsonic levels. During the super- 
sonic period of dissipation we expect that the cloud size 
will be decreasing - not so much due to stripping, but be- 
cause of its compaction. When the ram pressure, pv 2 , be- 
gins to decline, the cell will gradually transform from one 
contained by the shock shell to one contained by its own 
vorticity. This will be occassioned by some expansion, and 
as it expands, the cloud deceleration will increase, leading 
to a relatively rapid change of state. Though the density 
of the surrounding medium should by this time be rather 
large, it is likely — indeed it is imperative — that the den- 
sity in the central axis of the vortex ring should be orders 
of magnitude greater. The disappearance of the shock will 
allow the rapid cooling of cloud gas at the head of the 
shock, where it is densest. Yet the pressure on this gas 
is sustained by pressures from all directions: the gas near 
the axis is subjected to an axially symmetric centrifugal 
force from the vortex ring; toward the host galaxy there 
remains the ram pressure of the cloud at near sonic veloci- 
ties; and finally, pressure is exerted by material which has 
been stripped, but is re-joining its cloud, now falling with 
increased velocity due to the increased deceleration. 

To assess the plausibility and efficiency of star formation 
at this point we need to know the collapse time-scale for 
gas at the expected densities, the transit-time scale for gas 
in the axis of the vortex, and some notion of the duration 
of the vortex cell, which is responsible for maintaining the 
radial pressure along the axis of the cell. The time-scale 
for cloud collapse is, 

r ra!I «l/v^~4.7x 10 5 yr, (17) 

for the central density cited in §4.3. The time-scale for 
transit down an axis of 10 pc is ~ r e , given by Eq. 14, 
where, by assumption, v < c s . As the cell is slowed, this 
velocity will fall, in time providing a sure opportunity for 
gravitational collapse. Note that the cooling time-scale 
(Eq. 13) for baryon densities of order the central GC den- 
sities (n c i w 4 x 10 4 cm~ 3 ) is well within these time-scales. 
The cell itself should be stable for well over an eddy turn- 
over time-scale (Eq. 14), in keeping with the large dura- 
tions of more modest-sized vortex rings which may travel 
distances that are many times their diameters (Widnall & 
Sullivan 1973) before disruptionl. A more severe threat to 
the survival of the vortex cell is condensation of gas and 
star formation, but by then, we don't care. 

Star formation rates, and efficiencies are now of concern. 
In a study of 97 normal, and star-forming galaxies, Ken- 
nicutt (1977) found that the disk-averaged star formation 
rates were found to fit to a Schmidt (1959) law with index 
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N = 1.4 ± 0.15. It was found that, 

/ V \ 1-4±0.15 

S SFK = (2.5±0.7) x 10- 4 gas M Q yr- 1 kpc 

VlM pc 2 J 

(18) 

If we imagine gas with a density equal to peak central GC 
densities (§4.3), distributed over a 10 pc cube, then the 
indicative star formation density is ~ 0.25 M Q yr -1 pc~ 2 , 
and would yield 1.25 x 10 7 M within a plausible time- 
scale of 5 x 10 5 yr, a mass larger than the 10 6 M Q con- 
tained in the cube. Star formation efficiency would vary 
with index N —1, or, X°, a 4 . At projected final cloud densi- 
ties, this implies star formation efficiencies 40 times greater 
than that of a normal disk with a gas surface density of 
<~ 1 M Q pc~ 2 . This would rival that of central starbursts 
(Kennicutt 1977). 

5. IMPLICATIONS FOR FIELD GALAXY FORMATION 

A model in which galaxies are subjected to periodic 
accretion of low metallicity clouds is consistent with the 
"chaotic" scenario of galaxy formation of Scarlc & Zinn 
(1978), and Searle (1980), now often referred to as hi- 
erarchical structure formation. The GC formation sce- 
nario requires a picture in which galaxy growth is a pro- 
tracted process. That the accretors are probably DM- 
held, and generally relaxed, seems required by the spe- 
cific scenario developed here. Shull et al. (1996) estimate 
that these clouds have a total mass (baryons plus DM) 
of M to t — 10 9 ' 8 MqRiooT4. 3 , where R\qq is the radius 
in units of 100 kpc at which the column density of neu- 
tral hydrogen is 10 13 cm -2 , and T is in units of 10 4,3 K. 
This value is consistent with a mass distribution constant, 
fC f=s 10 5 M Q pc -1 , with a circular velocity v c ~ 20 kms" 1 , 
as observed by Blitz et al. (1998). Somewhat smaller 
clouds may also be realistic. 

The accretion of such a cloud would produce great 
turbulence. Might there be low-redshift counterparts of 
proto-late-type galaxies? High resolution VLA obser- 
vations of 5 actively star forming blue compact dwarfs 
(BCDs) reveal kinematically distinct clumps of H I, and 
turbulent outer H I envelopes (on scales of a fewx 100 pc). 
This is broadly suggestive of the accretion of H I clouds. 
That the median gas depletion time scale is on the order of 
1 Gyr also argues for extragalactic replenishment. Thus 
these young galaxies may resemble our Galaxy in its ear- 
liest formative stages. It is easily seen that the accretion 
of one or two hundred clouds of ~ 10 9 ' 8 M Q can account 
for the total mass of an L* galaxy within R gc w 100 kpc. 

While the formation of GCs is somewhat incidental to 
the process of galaxy formation as outlined above, yet 
they remain potentially sensitive probes of the state of 
the Galaxy at its early formative stages. The GC for- 
mation scenario suggests that among the requirements for 
cluster formation are two constraints on the (proto) galaxy 

- that it must have a large gas column density, and a 
mass distribution constant large enough to support re- 
quired circular velocities of J> 40 km s _1 , a value somewhat 
below that characteristic of dwarf spiral galaxies. The for- 
mer follows by the observation that ram pressure must 
slow the cloud to subsonic velocities: / p ga iv 2 S{t)dt w 

— / m c i(t) d ^^ dt, where p ga i is the galaxy baryon density, 
S(t) is the cloud surface area, and stripping will reduce the 



cloud baryonic mass, m c i, in time. Alternatively, 

». /;, s „,(«m«)s(«)« - -£ « ^ 

A detailed numerical simulation will be required to prop- 
erly evaluate these integrals, but clearly, if the column of 
galactic gas, J p(R)irr 2 l (R)dR is too small, then the equa- 
tion cannot be met. The mass-distribution requirement, 
together with the equation, V rot — VGK,, is essentially 
that K, £ 1.2 x H^Mgpc- 1 . This threshold for GC for- 
mation results from the observation that an infall velocity 
of 60 kms -1 , which is at the lower end of the range sat- 
isfying the energetic requirements for cloud compression, 
is ~ y/2 times the circular velocity at that distance. Once 
the galaxy mass is great enough, GCs would be formed. 
However, in response to an increased specific star forma- 
tion rate accompanying the formation of the bulge and 
the disk, the gas column should decrease, signaling an end 
of the halo and GC formation epoch. During the stage 
of this retreat, we would expect that, if it relaxed into a 
thick, slowly rotating disk, then GCs would have a lower 
probability of being formed if the cloud were plunging nor- 
mal to the disk, due to the lower gas columns. Thus, the 
metal-rich "disk" globulars, with their flattened distribu- 
tion (Armandroff 1993), and the thick disk itself, may be 
evidence favoring this interpretation. 

6. IMPLICATIONS FOR COSMOLOGY 

While the rcdshifts of formation implied by this sce- 
nario may strike the reader as quite low, recent findings, 
largely based on Hipparcos data, the inclusion of helium 
diffusion into the cores of stars, and an improved equation 
of state, yield new distance determinations based on the 
main sequence turnoff, which in turn imply that the oldest 
globular clusters may be much closer to 12 Gyr old. This is 
a ~ 20% reduction from pre- Hipparcos values (for a recent 
review, see Chaboyer 1998). If we allow that these oldest 
GCs were formed at z — 3, then given this GC age, we may 
derive H Q as a function of cosmology. Accordingly, we find 
that if Q = 1, we would require H a = 47.7 kms -1 Mpc -1 . 
If O = 0, H = 61.3 kms^Mpc" 1 . For a flat A model 
with A = 0.7, H = 67.8 kms -1 Mpc -1 . If the formation 
epoch instead were at z ~ 4, then the attributed Hub- 
ble constant is increased by less than 6.5% for all three 
models. 

It has been reported that GCs may have a significant 
range of ages. Recent estimates are in the range of ~ 5 Gyr 
(Fusi-Pecci et al. 1995; Chaboyer, Demarque & Saraje- 
dini 1996). Using standard astrophysical formulae, I find 
that the range of time between the plausible redshift range 
z = 1.5 — 3, the expected range of GC formation redshifts, 
is, 1.79, 2.44 and 2.20 Gyr respectively for the O = 1, 0, 
and A = 0.7 cosmologies, where for self-consistency I have 
used the specific value of H Q derived for each case. These 
ranges are smaller than the observed value cited above, 
but not so much as to cause alarm at this early stage. 

7. CONCLUSION 

High resolution HST spectra at low redshift have dis- 
closed the existence of a surprisingly large population of 
Lyman a clouds, most of which are thought to be clustered 
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within 2 /i^ 5 Mpc of bright field galaxies. The juxtaposi- 
tion of evidence for cloud clustering, and for major accre- 
tion events of low metallicity gas onto large field galaxies, 
suggests a causal relationship. By following the common- 
alities of cloud clustering among the low-redshift popula- 
tion out to high-redshift, it becomes clear that clouds at 
high-redshift certainly had the opportunity to accrete to 
the protogalaxies about which they are thought to clus- 
ter. The dramatic disappearance of Lyman a clouds at 
high-redshifts provides the link between the clustering of 
clouds and the heightened star formation rates in field 
galaxies seen at z ~ 1 — 3 (Madau et al. 1998). It is 
suggested that the central regions of many of these clouds 
might plausibly have been transformed into GCs, while 
the less-strongly held gas may have contributed to the 
formation of the stellar halo. It has been shown that the 



energetics of cloud compression are favorable, as are the 
numerical coincidences of the comoving densities of GCs, 
and clouds at high-redshift. The juxtaposition of the pro- 
jected redshifts of formation with new globular cluster ages 
implies values of H a that are reasonable. However, the pre- 
dicted cosmology-specific age range of GCs appears to be 
lower than recent work would imply. This interesting fact 
promises to be a goad to stimulate future work. 

I would like to thank the anonymous referee for many 
comments and suggestions that have improved this paper. 
In addition, I thank Hyron Spinrad, Daniel Stern, Christo- 
pher F. McKee, Robert Fisher, and Dean McLaughlin for 
helpful suggestions during the planning and production of 
this paper. 
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